Strong Desynchronizing Effects of Weak Noise 
in Globally Coupled Systems 
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In assemblies of globally coupled dynamical units, weak noise perturbing independently the individ- 
ual units can cause anomalous dispersion in the synchronized cloud of the units in the phase space. 
When the noise-free dynamics of the synchronized assembly is nonperiodic, various moments of the 
linear dimension of the cloud as a function of the noise strength exhibit multiscaling properties 
with parameter-dependent scaling exponents. Some numerical evidence of this peculiar behavior as 
well as its interpretation in terms of a multiplicative stochastic process with small additive noise is 
provided. Universality of the phenomenon is also discussed. 



I. INTRODUCTION 

Oscillatory behavior, periodic or nonperiodic, observ- 
able on a macroscopic scale may often be a result 
of collective synchronization of a population of micro- 
oscillators Consider, for simplicity, a population 
of identical oscillators with global coupling. Systems of 
globally and uniformally coupled oscillators in fact serve 
as a natural model for Josephson junction arrays 
multimode lasers |||7]] and various biological systems such 
as flashing swarms of fireflies and pacemaker cells 
responsible for circadian rhythms |H|. Under suitable 
conditions, the whole population will exhibit complete 
synchrony when the oscillators are free from intrinsic or 
extrinsic randomness. 

In practical situations, however, the oscillators would 
be more or less noisy or non-identical, so that, instead 
of exhibiting perfect synchrony, they will be distributed 
in the phase space to form a cloud of a finite extension. 
It seems quite natural to expect that, as far as the noise 
remains sufficiently weak, the long-time average of the 
linear dimension of this cloud, denoted by (r), should be 
proportional to the noise strength. This is actually the 
case when the oscillators are of the limit cycle type. We 
shall find that, when the oscillators are chaotic, in con- 
trast, such linear dependence breaks down under broad 
conditions, and is replaced by a more general power-law 
dependence. The main goal of the present paper is to 
explain theoretically why such peculiar behavior is pos- 
sible for globally coupled chaotic oscillators. This will be 
verified with some numerical evidence for a number of 
population models. Our theory, which is quite similar to 
the one developed earlier to explain a certain unique fea- 
ture of turbulence in nonlocally coupled systems [ ll|-^5| , 
suggests that the general moments (r q ) as a function of 
the noise strength exhibits a simple multiscaling law. As 
a natural consequence from the same theory, similar mul- 
tiscaling behavior is expected to occur even for the popu- 
lations of limit cycle oscillators provided that, apart from 
the random noise, the oscillators are subjected to another 
random force which is common over the whole popula- 
tion. The crucial point here is that the average motion 



has to be stochastic for the multiscaling behavior to oc- 
cur. 

In Sec.0, we start with a general class of differential- 
equation models for globally coupled elements with ad- 
ditive noise. Then, as a specific model belonging to this 
class, we study numerically a population of the chaotic 
Rossler oscillators , and show the anomaly in the size 
of the synchronized cloud of the oscillators in the phase 
space. In Sec. Ill, we develop a theory to explain the ori- 
gin of a general scaling law including the results of Sec.|| 
as a special case. Our theory becomes almost identi- 
cal with the one developed in Refs. [|l2| if the effects of 
the external noise are replaced by the effects of spatial 
nonuniformity. The arguments in Sec. Ill suggest that 
the phenomenon of concern should be so universal that 
a number of restrictions imposed on our model could be 
removed. We provide some evidence for this in Sec.[y|, 
where three different types of population model will be 
discussed. The final section summarizes our main con- 
clusions. 



II. MODEL AND NUMERICAL SIMULATION 

Consider a population of N identical units with the 
intrinsic dynamics given by X = F(X), where N is an 
arbitrary number. Typically, the dynamical units con- 
sidered are chaotic oscillators. Introducing all-to-all type 
linear coupling and also external additive noise driving 
the oscillators individually, we have the following system 
of coupled differential equations: 



dXj 
dt 



F(X t ) + K-[X(t)-X t (t)] 
(i = l,...,N). 



■f'Vi (*) 



(2.1) 



Here X is the simple average of Xi over the population, 
i.e., 



(2.2) 



»=i 



and if is a positive coupling constant, thus working in 
favor of synchronous behavior of the population; rj 4 (t) 



1 



represents the external noise with suitably normalized in- 
tensity applied independently on the oscillators, so that 
the coefficient / measures the intensity of the noise. For 
the sake of simplicity, the coupling term has been so ar- 
ranged that the dynamics of a given unit which is in per- 
fect synchrony with the average motion would be iden- 
tical with its own dynamics without coupling. The cou- 
pling strength K may be generalized to a coupling ma- 
trix, without qualitatively affecting the whole arguments 
which follow. 

In the absence of noise, the population will be attracted 
to a perfectly synchronized motion for sufficiently large 
K. In the one-oscillator phase space, the dynamics could 
then be imagined as an evolution of a single point. When 
a small random noise is introduced, the resulting dynam- 
ics can no longer be represented by a point in the phase 
space. Instead, we have a cloud of a finite extension 
whose shape will be changing variously in time. The 
effects of noise on the collective dynamics could conve- 
niently be characterized by the average linear dimension 
of the oscillator cloud in the phase space defined by 

1 N 

r(i) = -£|*(*)-*< (*)|- (2.3) 

i=l 

We define its various moments by a long-time average of 
(r»>: 

(r«)= Urn ± / r«(t)dt. (2.4) 

Our main concern below is the dependence of (r q ) on 
/, K and q. Our intuition may tell that (r q ) will sim- 
ply be proportional to f q , which seemed to be true from 
some numerical simulation jTi| . However, we see below 
that this is not always valid. Instead, (r q ) behaves under 
broad conditions like 

(r q ) oc f a{q) (2.5) 

for sufficiently small /, where the exponent a(q) is not 
linear in q and also changes continuously with K and 
other system parameters. 

The above power-law dependence can be confirmed by 
a numerical analysis of globally coupled Rossler oscilla- 
tors in the chaotic regime. By assuming the coupling 
matrix to be diagonal, the system is governed by the 
equations 

Xi = -(Yi + Zi) + K_(X -Xt)+f- m,x, 
Yi=X i + aY i +K ■ (Y - Y^+ f ■ i^ Y , 
Z t = b- cZi + XiZi +K-(Z-Zi)+f- m,z, 

(2.6) 

where the noise components r)i lV (y = X, Y, Z) are as- 
sumed to be mutually independent and white Gaussian 
with vanishing mean and unit variance. We will choose 
the standard parameter values a = 0.3, b = 0.2 and 



c = 5.7. It turns out that in the absence of noise the 
state of perfect synchrony is stable above some critical 
value K c of K, where K c ~ 0.4. We will work with this 
condition throughout, and concentrate on the the size of 
the oscillator cloud when the noise sets in. Some numer- 
ical results for the first moment (r) are summarized in 
Figs. and |. 
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FIG. 1. log (r) vs. log / calculated numerically from Eq. 
( ^.fi| ) with TV = 64 for several values of K. The top line 
corresponds to K — 0.142 and the bottom one to K = 0.170, 
with the uniform interval of 0.004. 




0.14 0.15 0.16 0.17 0.18 0.19 



coupling strength K 
FIG. 2. Exponent a(l) vs. coupling strength K, obtained 
from the data of Fig. |l] using the method of least-squares. 

In Figure [j], log-log plots of (r) versus / are displayed 
for different values of K. The power-law behavior as 
given by Eq. (|]^) with q = 1 seems to hold well with 
the exponent a(l) changing with K. Fig. || shows the 
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dependence of the estimated values of a(l) on K. For suf- 
ficiently strong coupling, the exponent saturates to the 
normal value a(l) = 1, while for weaker coupling a(l) 
varies between and 1. Since the noise is sufficiently 
small, the value of a(l) less than 1 implies anomalous 
amplification in the size of the oscillator cloud as com- 
pared with the normal case. 



III. ORIGIN OF ANOMALOUS FLUCTUATIONS 



In this section, we present a theory on the origin 
of anomalous size fluctuations of a synchronized cluster 
found numerically in the preceding section. Our theory 
can be developed quite in parallel with the theory devel- 
oped before on the multiscaled turbulence in nonlo call y 
coupled systems [[n] ^5|. Let us start with Eq. (2.1). 
Note that the average motion of the population obeys 
the equation 



dX 
~~dt 



FX) 



O 



N ^ 



(3.1) 



Xi (t) - X (t) is 



Thus, the deviation defined by a;, (t) 
governed by 

= (DF (X (*)) - K) Xi + f ■ m (t) + OQxtf), 

(3.2) 

where DF (X(t)) is the Jacobian of F(X) at X = 

X (t). For sufficiently large t, the vector Xi (t) is ex- 
pected to become parallel with the Lyapunov eigenvector 
of DF{X(t)) — K corresponding to the largest Lyapunov 
eigenvalue. By denoting the amplitude of this eigencom- 
ponent of Xi (t) as Xi (t), an d d ropping the suffix i for 
the sake of simplicity, Eq. ( |3.2[ ) is reduced to a scalar 
equation 



^ = X(t)x + f-r](t) + 0(x 2 ) , 



(3.3) 



where A(t) is the local Lyapunov exponent, and T)i(t) is 
the projection of rji(t) onto the corresponding Lyapunov 
eigen vect or. For the coupled Rossler oscillators given by 
Eq. (2.6), we have 



A(t) = X r (t)-K, 



(3.4) 



where A r (t) is the local Lyapunov exponent of the indi- 
vidual Rossler oscillator. Since the sign of the long-time 
average of A(i), which will be denoted by Ao, determines 
the stability of the perfect synchrony of the population in 
the absence of noise, the critical coupling strength for the 
coupled Rossler oscillators is simply given by K c = Ao. 

Since we are working with the condition that the long- 
time average of X(t) is negative, and also with sufficiently 
weak noise, there should exits a range of x where the 



nonlinear term of 0(x 2 ) as well as the noise term / • r](t) 
are negligible. Let this range be specified by x m i n <C x <C 
Xmaxi where x m i„ should be proportional to /. In this 
range, we have a linear multiplicative stochastic equation 



dx 
~dt 



\(t)x. 



(3.5) 



Note that the stochasticity here is of a deter min istic ori- 
gin. The stationary distribution from Eq. ( |3.5| ) can be 
discussed analytically in the following way. Introducing 
a new variable y by 



y = inM, 

we rewrite Eq. ([T^) as 
This is integrated to give 



y(t)-y(o) 



A (t) dt = A (t) 



(3.6) 



(3.7) 



(3.8) 



The time-dependent probability distribution function 
(PDF) of y, denoted by Q (y,t), satisfies 



Q{y,t) = / w(A,t)Q(y-A,0)dA, 



(3.9) 



where u> (A, t) is the PDF of A (t) and is equivalent to 
the transition probability between the states at times 
and t. From the definition of A, u> (A, t) is expected to 
converge to a Gaussian of mean t(Ao — K) (= t(K c — K)) 
and variance tD\ as t — » oo, where 



D A = lim - 




(A(t) - Xo)dt 



Thus, the stationary PDF of y becomes 
Q (y) oc e-^ , 

where 

2A 



(>0) . 



(3.10) 



(3.11) 



(3.12) 



Note that for the coupled Rossler oscillators, (3 is given 
by 



(3 = 



2(K - K c 
D x 



(3.13) 



The stationary PDF for the original variable x, denoted 
by P(x), must s atisfy the relation 2P(x)dx — Q(y)dy. 
Thus, from Eq. ( 3.11 ), we obtain 



P (x) OC X ^ ^ {Xmin < ^>~ X X ma x) 



(3.14) 
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Coming back to the original stochastic equation (3.3), 
the effects of additive noise and nonlinearity have now 
to be incorporated. There exists a characteristic value 
of x below which the additive noise term dominates the 
other terms. This value is what we denoted by x m i n - 
For x -C Xmin, the power-law divergence of P(x) will be 
saturated to a constant. On the other hand, there exists 
the second characteristic value of x, denoted by x max , 
above which the nonlinear term is the most dominant. If 
the nature of the nonlinearity is such that it decelerates 
rather than accelerates the growth of x, which we assume, 
the power-law decay of P(x) will be replaced by a much 
sharper decay above x max . As far as their dependence 
on / is concerned, z m >n an d x max may be specified as 



1. 



(3.15) 
(3.16) 



Thus, the entire profile of P(x) will be such that it 
is nearly constant below x = x m m, obeys the above- 
mentioned power law, and decays quickly above x = 
Xmax- On further idealization, we represent P(x) by the 
analytic form 



P(x) 




(0 < x < f) 
(f<x<l) 
(Kx) 



(3.17) 



where C is a normalization constant. 
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FIG. 3. Stationar y P DFs of P (f) calculated nu- 
merically from Eq. (2.6) for several values of noise 
strength /, where K — 0.16 is fixed. The curves 
correspond to / = 10 _1 (lowest level at lower £), 
lO- 1 - 25 ,^- 1 ' 5 ,^- 1 - 75 ,^- 2 ' 00 and 10~ 2 ' 25 (highest). 
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FIG. 4. Stationary PDFs P (£) calculated numerically from 
Eq. (2.6) for several values of coupling strength K, where 
/ = 0.01 is fixed. The curves correspond to K = 0.150(slowest 
decay), 0.155, 0.160, 0.165, 0.170, 0.175 and 0.185(fastest de- 
cay). 

Figures || and |i| show log- log plots of P(£) for the de- 
viation £ of the first component X obtain ed n umerically 
from the coupled Rossler oscillators (Eq. (51])) for some 
values of / and K . It is seen that each curve is approx- 
imately composed of three parts corresponding to the 
three characteristic regimes in Eq. (3.17). The range of 
£ where P(£) is nearly constant e xtend s proportionally 
to /, which is consistent with Eq. ( 3.17 ), while the sharp 
drop of P(£) seems to occur at som e £ in dependent of /, 
which is again consistent with Eq. ( |3. 17 ) with x replaced 

by e 

Up to this point, we have discussed the statistics of 
the deviation a; of a single element from the mean motion 
of the population. In what follows, we will identify the 
statistics of |x| with that of r, i.e., the deviation averaged 
over the whole population. As far as some qualitative 
features such as the power-law dependence of (r) on / is 
concerned, this assumption seems to be justified because 
the dynamical units are driven by a common multiplier 
X(t) by virtue of the g lobal nature of the coupling. 

By using Eq. ( 3.17 ) , it is thus straightforward to cal- 
culate various moments of r. As for the first moment, we 
have 



<r) 



\x\ P (x) dx oc 



f 1 



(0 < < 1) 
08 >1) 



(3.18) 



Comparing the above equation with the expression in 
Eq. (3.13) for the coupled Rossler oscillators, one may 
now understand the reason for the observed anomalous 
power-law dependence of (r) on the external noise when 
the system is not too far from the critical point K = K c . 

The observed change in the scaling exponent with K 
shown in Fig. @ still deviates considerably from Eq. 
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( |3 .18 ) with j3 given by Eq. (3.13). Specifically, the nu- 
merical results do not exhibit sharp changes near /3 = 
and f3 = 1. Such discrepancy seems to be d ue to the 
fact that the range of validity of Eq. ( 3.18 ) in terms 
of / shrinks to zero as (3 approaches or 1. A little 
more careful analysis shows that under fixed /, we have 
(r) - 1/| In /| as (3 -> 0, and (r) - |/"ln/| as /3 -> 1. 

Similar calculations for the higher moments (r 9 ) are 
also straightforward, and the results are simply given by 



A. Case of globally coupled maps 



(r q ) oc 



f" 



(0<P<q) 



(3.19) 



except for the weak logarithmic singularity mentioned 
above. Thus, the anomalous fluctuations could be vis- 
ible through higher moments in the range of stronger 
coupling where no anomaly is visible through lower mo- 
ments. Specifically, for the coupled Rossler oscillators, 
the range of K where the q-th moment behaves anoma- 
lously is given by 



K C <K<K C 



qD x 



(3.20) 



Near the critical point at w hich the average Lyapunov 
exponent vanishes, Eq. (3.3) is of the same form as the 
equation employed for discussing the so-called on-off in- 
termittency 18-21] with noise [ p2|j2^ ]. Since the dynam- 
ics of r(t) would qualitatively be the same as that of x(t) 
of a representative oscillator, this implies that noisy on- 
off intermittency could also be observed in r(t). This was 
confirmed numerically for the coupled Rossler oscillators, 
though we will not report its details. The only thing to 
be remarked is that the origins of the power law in the 
on-off intermittency and that of our present concern are 
completely independent. 



IV. UNIVERSALITY OF ANOMALOUS 
FLUCTUATIONS 

Up to the preceding section, our discussion has been 
based on the following three assumptions: 

1 . The system is described by a set of ordinary differ- 
ential equations 

2. The source of randomness working against com- 
plete synchrony is represented by external additive 
noise 

3. The constituents of the population are intrinsically 
chaotic 

In what follows, we will show that neither of these as- 
sumptions is necessary, which implies that the anomalous 
behavior of concern would be quite general. 



Discrete-time analogue to Eq. (2.1) is the system of 
globally coupled maps |24|,|25| with noise. Assuming for 
the sake of simplicity that the individual map is one- 
dimensional, we are concerned with the model equation 
of the form 

Xi (n + 1) = (1 - K) M(Xi (n)) + KM (xl^fj + fm 

(i = l,...,N). (4.1) 

One may develop arguments similar to the case of 
continuous-t ime dynamics, and derive a discretized ver- 
sion of Eq. (3.3): 



+ f- Vn + 0(xl). 



(4.2) 



Furthermore, the argument s lead ing to a stationary PDF 
similar to the form of Eq. (3.17) are almost the same as 
before. As an example, let us consider a generalized tent 
map for M(x) defined by 



M(x) 



x/a (0 < x < a) 

(l-x)/(l-a) (a<x<l) 



(4.3) 



The crucial condition for the occurrence of an oma ly of 
(r) is that the local Lyapunov exponent in Eq. (4.2) fluc- 
tuates between positive and negative values. Thus, the 
conventional tent map (a = 0.5) for which A = In 2 iden- 
tically is ruled out. I n a suitable range of a and K, the 
system given by Eq. (4.1) is confirmed to ex hibit power- 
law fluctuations in the form of Eq. ( |3.18 ). Estimated 
exponent a(l) of the first moment of r as a function of 
K is displayed in Fig. [|. 
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0.44 0.46 0.48 0.5 0.52 0.54 0.56 

coupling strength K 
FIG. 5. Exponent a (1) vs. coupling strength K calculated 
numerically from the generalized tent map (Eq. (4.3)) with 
a = 0.75 



B. Effects of heterogeneity 

Any kind of heterogeneity present in the population 
could give a source of incompleteness in synchrony. We 
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focus on the following inhomogeneous system of elements 
without additive noise. 



dXj 



= F{X l ,a l ) + K-(X (t)-Xi (t)) 



(4.4) 



Here the heterogeneity is represented by a weak distri- 
bution in a parameter a, its value for the i-th element 
being denoted by a^. The dynamics of the deviation x of 
a representative oscillator is described by 



^ = {DF(XU).n)-K)x 



da 



OF (X (t),a) 
da 



O \x\\5a- 



(4.5) 



where suffix i is dropped, and 5a denotes the deviation of 
a from its population average a. The effect of the second 
term on the right-hand side is esse ntially the same as that 
of the additive noise in Eq. (3.2), thus resulting in the 
anomalous power-law dependence of {r q ) on the strength 
of inhomogeneity. 



C. Populations of non-chaotic units 

We have seen that the anomalous fluctuations exhib- 
ited by the synchronized cluster can arise only if the lo- 
cal Lyapunov exponent associated with the individual 
unit fluctuates randomly. This means that we can never 
expect such anomaly for populations of non-chaotic dy- 
namical units. Even for this class of systems, however, 
the anomalous fluctuations may arise provided the whole 
population is driven externally by a common random 
force apart from the random noise considered previously. 
The system of this class will take the form 



dXj 
dt 



= F (X,) + K(X (t) - X t (*)) + G (t) + f Vi (t) , 

(4.6) 



where the individual dynamics X = F(X) is assumed 
to be non-chaotic, and G (t) represents the coherent ran- 
dom force independent of i. By virtue of the G(t) term, 
the local Lyapunov exponent now fluctuates randomly, 
so that the anomalous fluctuation in r could be recov- 
ered. This can be demonstrated with a population of 
oscillatory/excitable units of the FitzHugh-Nagumo type 
p6|p7| . The specific form of the model is given by 

^ = - e (X i -Xf-Y i )+K.(X-X i ) 

+ G sin <j)(t) + f ■ m,x{t), 

dY 

= aXi + b + K-(Y-Yj)+f- Vl , Y {t). (4.7) 

where Goshi(f>(t) represents external random force with 
4> generated from the dynamics of a random walker: 



dt 2 



(4. 



Some numerical results of the above model are shown in 
Figs. ^ and ^ where the parameter values are so cho- 
sen that the individual unit is non-oscillatory but ex- 
citable. Nontrivial power-law dependence of (r) on / 
with parameter-dependent exponent is again confirmed. 
Similar behavior of (r) of course persists when the indi- 
vidual dynamics becomes oscillatory. 
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logio(/) 

FIG. 6. log (r) vs. log / for several values of K, from 
0.01 (top) to 0.19 (bottom) with the unifor m i ncreme nt o f 
0.002, calculated numerically from Eq. (4.7) and (4.8); 
N = 64, (a,b,e) = (1.0,0.58,0.1), (Go, 7, g) = (0.5, 1.0,5.0). 
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coupling strength K 
FIG. 7. Exponent a(l) vs. coupling strength K obtained 
from the data of Fig.M using the method of least-squares. 
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V. SUMMARY 

In the present paper, we argued that in globally cou- 
pled systems of nonlinear dynamical units small noise 
or imposed heterogeneities can generally cause anoma- 
lously strong dispersion of the synchronized clusters in 
the phase space. This was demonstrated numerically 
with a number of population models with chaotic dy- 
namical units. The numerical results were explained the- 
oretically in terms of a multiplicative stochastic process 
with additive noise. It turned out that the crucial con- 
dition for the occurrence of such anomaly is the random 
fluctuations of the local Lyapunov exponent associated 
with the individual units. This fact suggested some pos- 
sible generalizations of the class of systems capable of 
exhibiting similar behavior. In particular, chaotic nature 
of the individual dynamics seemed unnecessary provided 
the population is subjected to a common random drive, 
and this was actually demonstrated with the population 
of the FitzHugh-Nagumo type excitable units. 
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